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Abstract 

We  propose  the  design  of  innovative  routing,  planning  and  coordination  strategies  for  robotic  networks 
and  their  application  to  Army  and  DoD  scenarios.  The  key  technological  challenge  is  the  decision  of  who 
does  what,  when  or,  equivalently,  how  are  tasks  partitioned  among  robots,  in  what  order  are  they  to  be 
performed,  and  along  which  routes  do  individual  robots  move. 

The  fundamental  novelties  and  our  recent  breakthroughs  motivating  this  proposal  are  manifold:  (1)  the 
application  of  queueing  theory  and  combinatorial  techniques  to  network  of  autonomous  robots  leads  to  a 
wide  range  of  new  relevant  problems,  (2)  novel  coordination  schemes  promise  to  successfully  achieve  various 
optimization  objectives  relying  only  upon  asynchronous  and  asymmetric  communication,  (3)  the  increasingly 
weaker  assumptions  imposed  on  routing  and  coordination  algorithms  are  rendering  them  practical  and  widely 
applicable. 

This  proposal  addresses  multi-dimensional  problems  of  relevance  in  Engineering  and  Computer  Science 
by  unifying  fundamental  concepts  from  multiple  domains  (robotics,  autonomy,  combinatorics,  and  network 
science).  The  proposal  aims  to  bridge  multiple  scientific  disciplines,  including  control  theory  and  theoretical 
computer  science  and  their  applications  to  multi-agent  systems,  robotics  and  sensor  networks. 

Throughout  the  duration  of  this  grant,  we  made  solid  contributions  in  numerous  key  areas: 

1.  stochastic  surveillance  for  quickest  detection, 

2.  dynamic  vehicle  routing, 

3.  coverage  control  and  equitable  partitioning, 

4.  knapsack  optimization  with  application  to  mixed  teams, 

5.  controllability  for  complex  networks, 

6.  cooperative  patrolling  for  robotic  and  camera  networks, 

7.  distributed  optimization  and  task  assignment 

In  the  following,  we  review  our  scientific  progress  and  accomplishments  in  each  of  these  areas. 


1 


1  Multi-robot  stochastic  surveillance  for  quickest  detection 

In  this  first  work  (see  (Jl)  below),  we  design  persistent  surveillance  strategies  for  the  quickest  detection  of 
anomalies  taking  place  in  an  environment  of  interest.  From  a  set  of  predefined  regions  in  the  environment,  a 
team  of  autonomous  vehicles  collects  noisy  observations,  which  a  control  center  processes.  The  overall  objective 
is  to  minimize  detection  delay  while  maintaining  the  false  alarm  rate  below  a  desired  threshold.  We  present  joint 
(i)  anomaly  detection  algorithms  for  the  control  center  and  (ii)  vehicle  routing  policies.  For  the  control  center, 
we  propose  parallel  cumulative  sum  (CUSUM)  algorithms  (one  for  each  region)  to  detect  anomalies  from  noisy 
observations.  For  the  vehicles,  we  propose  a  stochastic  routing  policy,  in  which  the  regions  to  be  visited  are 
chosen  according  to  a  probability  vector.  We  study  stationary  routing  policy  (the  probability  vector  is  constant) 
as  well  as  adaptive  routing  policies  (the  probability  vector  varies  in  time  as  a  function  of  the  likelihood  of  regional 
anomalies).  In  the  context  of  stationary  policies,  we  design  a  performance  metric  and  minimize  it  to  design  an 
efficient  stationary  routing  policy.  Our  adaptive  policy  improves  upon  the  stationary  counterpart  by  adaptively 
increasing  the  selection  probability  of  regions  with  high  likelihood  of  anomaly.  Finally,  we  show  the  effectiveness 
of  the  proposed  algorithms  through  numerical  simulations  and  a  persistent  surveillance  experiment. 

In  our  more-recent  work  (see  (J2)  below),  we  provide  analysis  and  optimization  results  for  the  mean  first 
passage  time ,  also  known  as  the  Kemeny  constant ,  of  a  Markov  chain.  First,  we  generalize  the  notion  of  Kemeny 
constant  to  environments  with  heterogeneous  travel  and  service  times,  denote  this  generalization  as  the  weighted 
Kemeny  constant ,  and  we  characterize  its  properties.  Second,  for  reversible  Markov  chains,  we  show  that  the 
minimization  of  the  Kemeny  constant  and  its  weighted  counterpart  can  be  formulated  as  convex  optimization 
problems  and,  moreover,  as  semidefinite  programs.  Third,  we  apply  these  results  to  the  design  of  stochastic 
surveillance  strategies  for  quickest  detection  of  anomalies  in  network  environments.  We  numerically  illustrate 
the  proposed  design:  compared  with  other  well-known  Markov  chains,  the  performance  of  our  Kemeny-based 
strategies  are  always  better  and  in  many  cases  substantially  so.  These  results  leave  open  numerous  directions 
for  further  research.  First,  we  designed  surveillance  policy  only  for  single  agent  systems  and  it  would  be  of 
practical  interest  to  consider  the  case  where  there  are  multiple  agents.  Second,  it  would  be  useful  to  understand 
bounds  on  the  design  of  of  the  mean  first  passage  time  for  general  graph  topologies.  Finally,  we  treat  only  the 
optimization  of  the  transition  matrix  of  the  graph.  It  would  be  of  interest  to  study  how  we  can  optimize  the 
weight  matrix  W  in  conjunction  with  the  transition  matrix.  This  can  have  the  interpretation  of  optimizing  the 
’’capacity”  or  ’’resistance”  of  the  graph,  a  topic  in  optimization  which  is  of  independent  interest. 

Finally,  in  (J3),  we  provide  generalized  notions  and  analysis  methods  for  the  hitting  time  of  random  walks  on 
graphs.  The  hitting  time,  also  known  as  the  Kemeny  constant  or  the  mean  first  passage  time,  of  a  random  walk 
is  widely  studied,  however,  only  limited  work  is  available  for  the  multiple  random  walker  scenario.  In  this  work 
we  provide  a  novel  method  for  calculating  the  hitting  time  for  a  single  random  walker  as  well  as  the  first  analytic 
expression  for  calculating  the  hitting  time  for  multiple  random  walkers,  which  we  denote  as  the  group  hitting 
time.  We  also  provide  closed  form  solution  for  calculating  the  hitting  time  between  specified  nodes  for  both 
the  single  and  multiple  random  walker  cases.  Our  results  allow  for  the  multiple  random  walks  to  be  different 
and,  moreover,  for  the  random  walks  to  operate  on  different  subgraphs.  Finally,  using  sequential  quadratic 
programming,  we  show  the  combination  of  transition  matrices  that  generate  the  minimal  group  hitting  time  for 
various  graph  topologies  are  often  different. 
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2  Dynamic  vehicle  routing 

In  this  thrust,  we  completed  the  survey  submission  (Jl)  below.  This  paper  surveys  recent  concepts  and  algo¬ 
rithms  for  dynamic  vehicle  routing  (DVR),  that  is,  for  the  automatic  planning  of  optimal  multi- vehicle  routes 
to  perform  tasks  that  are  generated  over  time  by  an  exogenous  process.  We  consider  a  rich  variety  of  scenarios 
relevant  for  robotic  applications.  We  begin  by  reviewing  the  basic  DVR  problem:  demands  for  service  arrive  at 
random  locations  at  random  times  and  a  vehicle  travels  to  provide  on-site  service  while  minimizing  the  expected 
wait  time  of  the  demands.  Next,  we  treat  different  multi-vehicle  scenarios  based  on  different  models  for  de¬ 
mands  (e.g.,  demands  with  different  priority  levels  and  impatient  demands),  vehicles  (e.g.,  motion  constraints, 
communication  and  sensing  capabilities),  and  tasks.  The  performance  criterion  used  in  these  scenarios  is  either 
the  expected  wait  time  of  the  demands  or  the  fraction  of  demands  serviced  successfully.  In  each  specific  DVR 
scenario,  we  adopt  a  rigorous  technical  approach  that  relies  upon  methods  from  queueing  theory,  combinatorial 
optimization  and  stochastic  geometry.  First,  we  establish  fundamental  limits  on  the  achievable  performance,  in¬ 
cluding  limits  on  stability  and  quality  of  service.  Second,  we  design  algorithms,  and  provide  provable  guarantees 
on  their  performance  with  respect  to  the  fundamental  limits. 

In  the  article  (J2),  we  considered  the  problem  of  dynamic  vehicle  routing  under  exact  time  constraints  on 
servicing  demands.  Demands  for  service  are  generated  in  an  environment  as  follows:  uniformly  randomly  in 
space  and  Poisson  in  time.  Every  demand  needs  to  be  serviced  exactly  after  a  fixed,  finite  interval  of  time  after 
it  is  generated.  We  design  routing  policies  for  a  service  vehicle  to  maximize  the  fraction  of  demands  serviced 
at  steady-state.  The  main  contributions  are  as  follows.  First,  we  demonstrate  that  this  problem  is  described 
by  an  appropriate  directed  acyclic  graph  structure  which  leads  to  a  computationally-efficient  routing  algorithm 
based  on  a  longest-path  computation.  Second,  we  provide  two  analytic  lower  bounds  on  the  service  fraction 
of  the  longest  path  policy.  The  first  bound  is  relative  to  an  optimal,  non-causal  version  of  the  policy,  i.e. ,  a 
policy  based  on  knowledge  of  all  future  demand  requests.  The  second  bound  is  an  explicit  function  of  demand 
generation  rate,  and  therefore,  useful  as  a  design  tool.  We  also  present  numerical  results  to  support  the  analytic 
bounds  and  to  shed  light  on  parameter  regimes  where  the  analytic  bounds  are  not  conclusive. 

In  (J3),  we  address  a  class  of  heterogeneous  multi-vehicle  task  assignment  and  routing  problem.  We  pro¬ 
pose  two  distributed  algorithms  based  on  gossip  communication:  the  first  algorithm  is  based  on  a  local  exact 
optimization  and  the  second  is  based  on  a  greedy  heuristic.  We  consider  the  case  where  a  set  of  heterogeneous 
tasks  arbitrarily  distributed  in  a  plane  has  to  be  serviced  by  a  set  of  robots,  each  with  a  given  movement  speed 
and  task  execution  speed.  Our  goal  is  to  minimize  the  maximum  execution  time.  We  propose  two  distributed 
and  asynchronous  algorithms:  the  first  one  is  based  on  the  iterative  optimization  of  the  local  task  assignment 
between  pairs  of  vehicles,  the  second  one  is  based  on  local  task  exchange  of  assigned  tasks,  one  by  one,  be¬ 
tween  couples  of  vehicles.  For  both  algorithms  we  provide  deterministic  bounds  to  their  performance.  The 
proposed  approaches  are  distributed  algorithms  easy  to  implement  in  a  networked  system  and  have  favorable 
computational  complexity  with  respect  to  the  ratio  k/n  between  the  number  of  tasks  and  vehicles. 
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3  Coverage  control  and  equitable  partitioning 

In  applications  such  as  environmental  monitoring  or  warehouse  logistics  a  team  of  robots  is  asked  to  perform 
tasks  over  a  large  space.  The  distributed  environment  partitioning  problem  consists  of  designing  control  and 
communication  laws  for  individual  robots  such  that  the  team  divides  a  space  into  regions  in  order  to  optimize  the 
quality  of  service  provided.  Coverage  control  additionally  optimizes  the  positioning  of  robots  inside  of  a  region. 
We  consider  a  comprehensive  set  of  scenarios  with  different  communication  protocols  and  spatial  constraints. 

In  recent  work  (see  Cl  and  J1  below),  we  consider  the  problem  of  optimal  coverage  with  area-constraints 
in  a  mobile  multi-agent  system.  For  a  planar  environment  with  an  associated  density  function,  this  problem 
is  equivalent  to  dividing  the  environment  into  optimal  subregions  such  that  each  agent  is  responsible  for  the 
coverage  of  its  own  region.  For  this  work,  we  design  a  continuous-time  distributed  policy  which  allows  a  team  of 
agents  to  achieve  a  convex  area-constrained  partition  of  a  convex  workspace.  Our  work  is  related  to  the  classic 
Lloyd  algorithm,  and  makes  use  of  generalized  Voronoi  diagrams.  We  also  discuss  practical  implementation  for 
real  mobile  networks. 
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4  Knapsack  optimization  with  application  to  mixed  teams 

In  a  recent  journal  article  (Jl)  below,  we  study  a  class  of  non-convex  optimization  problems  involving  sigmoid 
functions.  Examples  of  sigmoid  utility  functions  include  the  correctness  of  human  decisions  as  a  function  of 
the  decision  time,  the  effectiveness  of  human-machine  communication  as  a  function  of  the  communication  rate, 
human  performance  in  multiple  target  search  as  a  function  of  the  search  time,  advertising  response  as  a  function 
of  the  investment,  and  the  expected  profit  in  bidding  as  a  function  of  the  bidding  amount. 

We  show  that  sigmoid  functions  impart  a  combinatorial  element  to  the  optimization  variables  and  make 
the  global  optimization  computationally  hard.  We  formulate  versions  of  the  knapsack  problem,  the  generalized 
assignment  problem  and  the  bin-packing  problem  with  sigmoid  utilities.  We  merge  approximation  algorithms 
from  discrete  optimization  with  algorithms  from  continuous  optimization  to  develop  approximation  algorithms 
for  these  NP-hard  problems  with  sigmoid  utilities. 
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5  Controllability  for  complex  networks 

In  a  recent  technical  report  we  studied  the  problem  of  controlling  complex  networks,  that  is,  the  joint  problem 
of  selecting  a  set  of  control  nodes  and  of  designing  a  control  input  to  steer  the  network  to  a  target  state.  For 
this  problem  (i)  we  propose  a  metric  to  quantify  the  difficulty  of  the  control  problem  as  a  function  of  the 
required  control  energy,  (ii)  we  derive  bounds  based  on  the  system  dynamics  (network  topology  and  weights) 
to  characterize  the  tradeoff  between  the  control  energy  and  the  number  of  control  nodes,  and  (iii)  we  propose 
a  distributed  strategy  with  performance  guarantees  for  the  control  of  complex  networks.  In  our  strategy  we 
select  control  nodes  by  relying  on  network  partitioning,  and  we  design  the  control  input  by  leveraging  optimal 
and  distributed  control  techniques.  Our  findings  show  for  instance  that  (i)  if  the  number  of  control  nodes  is 
constant,  then  the  control  energy  increases  exponentially  with  the  number  of  the  network  nodes,  (ii)  if  the 
number  of  control  nodes  is  a  fixed  fraction  of  the  network  nodes,  then  certain  networks  can  be  controlled  with 
constant  energy  independently  of  the  network  dimension,  and  (iii)  clustered  networks  may  be  easier  to  control 
because,  for  sufficiently  many  control  nodes,  the  control  energy  depends  only  on  the  controllability  properties  of 
the  clusters  and  on  their  coupling  strength.  We  validate  our  results  with  examples  from  power  networks,  social 
networks,  and  epidemics  spreading. 
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6  Cooperative  patrolling  for  robotic  and  camera  networks 

Coordinated  teams  of  autonomous  agents  have  recently  been  used  for  many  tasks  requiring  repetitive  execution, 
including  the  monitoring  of  oil  spills,  the  detection  of  forest  fires,  the  track  of  border  changes,  and  the  patrol 
(surveillance)  of  an  environment.  The  surveillance  of  an  area  of  interest  requires  the  agents  to  continuously  and 
repeatedly  sweep  the  environment,  and  the  challenging  problem  consists  in  scheduling  the  agents  trajectories 
so  as  to  optimize  a  certain  performance  criteria. 

Our  research  in  this  direction  has  considered  the  patrolling  of  an  area  of  interest  by  means  of  a  team  of  mobile 
robots  and  by  means  of  a  network  of  autonomous  cameras.  For  both  cases,  we  define  suitable  performance 
functions,  we  characterize  optimal  trajectories,  and  we  design  centralized  and  distributed  algorithms  to  steer 
the  agents  along  optimal  trajectories. 

Regarding  the  patrolling  problem  with  mobile  robots,  in  article  (Jl)  below,  we  assume  the  robots  to  be 
identical  and  capable  of  sensing  and  communicating  within  a  certain  spatial  range,  and  of  moving  with  bounded 
speed.  We  represent  the  environment  as  a  graph,  in  which  the  vertices  correspond  to  physical  and  strategically 
important  locations,  and  in  which  the  edges  denote  the  possibility  of  moving  and  communicating  between 
locations.  Regarding  the  performance  criteria  of  a  patrolling  trajectory,  we  consider  (i)  the  time  gap  between 
any  two  visits  of  the  same  region,  called  refresh  time,  and  (ii)  the  time  needed  to  inform  the  team  of  robots 
about  an  event  occurred  in  the  environment,  called  latency.  Our  contributions  to  the  patrolling  problem  with 
mobile  robots  are  as  follows.  First,  we  mathematically  formalize  our  patrolling  problem  with  refresh  time  and 
latency  cost  functions.  We  show  that,  in  general,  designing  team  trajectories  with  minimum  refresh  time  is 
an  NP-  hard  optimization  problem,  and  we  identify  cases  (for  instance  when  the  roadmap  representing  the 
environment  has  a  chain  or  tree  structure)  in  which  an  optimal  trajectory  can  be  computed  efficiently.  Second, 
we  propose  different  approximate  trajectories  for  general  roadmaps,  and  we  characterize  the  performance  of 
our  approximate  patrolling  trajectories.  Third,  we  develop  distributed  algorithms  for  the  robots  to  synchronize 
along  the  proposed  approximate  trajectories.  Our  algorithms  leverage  on  tools  from  convex  and  combinatorial 
optimization,  and  rely  upon  mild  robots  communication  assumptions.  Fourth  and  finally,  we  validate  our  findings 
through  simulations  and  experiments:  we  use  the  Player/Stage  simulation  software  to  show  the  effectiveness 
and  the  robustness  of  our  second  patrolling  procedure  in  a  campus  environment,  and  we  conduct  an  experiment 
with  real  hardware  in  an  indoor  environment  with  obstacles.  The  experiments  confirm  the  robustness  of  the 
proposed  strategies  against  noise  and  unmodeled  dynamics. 

Regarding  the  patrolling  problem  with  a  network  of  autonomous  cameras,  in  article  (J2)  below,  we  consider 
Pan-Tilt-Zoom  (PTZ)  cameras  installed  at  important  locations.  We  assume  the  cameras  to  move  their  field 
of  view  (f.o.v.)  to  cooperatively  surveil  the  whole  environment.  As  for  the  case  of  robots,  we  model  the 
environment  with  a  robotic  roadmap.  We  develop  algorithms  for  the  cameras  to  self-organize  and  to  detect 
intruders  in  the  environment,  that  appear  at  arbitrary  locations  and  times.  We  consider  static  intruders,  which 
remain  stationary,  and  dynamic  intruders,  which  move  to  avoid  detection,  if  possible.  As  performance  criteria 
we  consider  the  worst-case  and  the  average  detection  times,  that  is  the  longest  time  and  the  average  time 
needed  for  the  cameras  to  detect  intruders.  Our  contributions  in  (J2)  are  as  follows.  First,  we  define  the  camera 
surveillance  problem  for  the  detection  of  static  and  dynamic  intruders,  and  we  formalize  the  worst-case  and 
average  detection  times  of  static  and  dynamic  intruders.  Second,  we  exhaustively  discuss  the  case  of  static 
intruders.  We  show  that,  for  tree  and  ring  roadmaps,  cameras  trajectories  with  minimum  worst-case  detection 
time  can  be  designed  by  solving  a  continuous  graph  partitioning  problem.  For  general  cyclic  roadmaps,  our 
trajectories  based  on  continuous  partitions  are  proved  to  be  optimal  up  to  a  factor  2.  Third,  for  the  case  of 
dynamic  intruders,  we  derive  a  necessary  and  sufficient  condition  on  the  cameras  locations  for  the  existence 
of  a  trajectory  with  finite  detection  time.  We  focus  on  ring  and  tree  roadmaps.  In  particular,  for  the  case 
of  ring  roadmaps  we  design  a  trajectory  with  detection  time  within  a  factor  3/2  of  optimal.  Instead,  for  tree 
roadmaps,  the  performance  of  our  trajectory  is  within  a  factor  2  of  optimal.  Fourth  and  finally,  we  consider 
three  different  communication  models,  and  we  propose  distributed  algorithms  for  the  cameras  for  continuous 
graph  partitioning  in  all  these  scenarios.  In  particular:  our  first  algorithm  assumes  a  synchronous  mode  of 
operation  of  the  cameras;  our  second  algorithm  assumes  an  asymmetric  broadcast  communication  model  and 
extends  the  class  of  block-coordinate  descent  algorithms  to  the  constrained  case;  and  our  third  algorithm  only 
requires  gossip  communication.  We  prove  convergence  of  all  these  algorithms,  and  we  analyze  their  performance 
in  a  simulation  study. 
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7  Distributed  optimization  and  task  assignment 

The  increasing  interest  in  performing  complex  tasks  via  multi-agent  systems  (e.g.,  sensor  and  robotic  networks) 
has  raised  the  interest  in  solving  distributed  optimization  problems.  The  fundamental  paradigms  in  distributed 
computation  are  that:  (i)  information,  relevant  for  the  solution  of  the  problem,  is  distributed  all  over  a  network 
of  processors  with  limited  memory  and  computation  capability,  and  (ii)  the  overall  computation  relies  only  on 
local  computation  and  information  exchange  amongst  neighboring  processors. 

Our  research  in  this  direction  has  focused  on  two  complementary  lines.  In  the  first  line  of  research  (see 
journal  article  (Jl)  below),  we  consider  a  distributed  version  of  linear  programs.  Optimizing  linear  objectives 
over  linear  constraints  takes  a  central  role  in  the  optimization  literature  and  thus  deserves  particular  attention 
also  in  distributed  computation.  Each  processor  in  the  network  is  assigned  only  the  information  relative  to  a 
subset  of  the  decision  variables.  The  objective  is  to  reach  an  agreement  on  a  global  minimum  of  the  problem, 
if  one  exists,  or  to  agree  that  the  problem  is  either  unbounded  or  infeasible.  We  show  in  this  paper  that  the 
Constraints  Consensus  algorithm  is  for  certain  problems  dual  to  our  new  algorithm.  For  the  new  algorithm,  we 
consider  linear  programs  in  equality  form,  which  are  allowed  to  be  fully  degenerate.  The  Constraints  Consensus 
can  in  general  not  be  applied  to  such  problems.  With  respect  to  this  we  propose  distributed  assignment  problems 
as  an  important  problem  class,  which  can  be  efficiently  solved  by  our  new  algorithm.  The  new  formulation  allows 
us  to  adopt  all  the  tools  known  in  the  classical  simplex  literature  to  this  problem  setup.  In  particular  simple  and 
efficient  decision  rules  for  the  pivot  iteration  can  be  provided  which  are  a  simple  mean  to  deal  with  degenerate 
linear  programs.  Finally,  our  new  algorithm  can  deal  with  infeasible  and  unbounded  problems. 

In  the  second  line  of  research,  see  journal  article  (J2)  below,  we  consider  distributed  approaches  for  the 
solution  of  random  convex  programs  (RCP).  RCPs  are  convex  optimization  problems  with  a  (usually  large) 
number  N  of  randomly  extracted  constraints;  they  arise  in  several  application  areas,  especially  in  the  context 
of  decision  under  uncertainty.  We  here  consider  a  setup  in  which  instances  of  the  random  constraints  (the 
scenario)  are  not  held  by  a  single  centralized  processing  unit,  but  are  instead  distributed  among  different  nodes 
of  a  network.  Each  node  ’’sees”  only  a  small  subset  of  the  constraints,  and  may  communicate  with  neighbors. 
The  objective  is  to  make  all  nodes  converge  to  the  same  solution  as  the  centralized  RCP  problem.  To  this 
end,  we  develop  two  distributed  algorithms  that  are  variants  of  the  constraints  consensus  algorithm:  the  active 
constraints  consensus  (ACC)  algorithm,  and  the  vertex  constraints  consensus  (VCC)  algorithm. 

The  major  contributions  of  this  work  are  threefold.  First,  we  develop  two  novel  algorithms,  ACC  and  VCC, 
for  the  computation  of  distributed  solution  to  RCP.  We  show  that  the  ACC  algorithm  converges  to  the  global 
solution  in  finite  time,  and  that  it  requires  almost  surely  bounded  communication  at  each  iteration.  We  give  some 
numerical  evidence  of  the  fact  that  the  ACC  algorithm  converges  in  a  number  of  iterations  that  is  linear  in  the 
communication  graph  diameter.  We  also  provide  numerical  evidence  that  parallel  implementation  of  the  ACC 
algorithm  significantly  reduces  the  computation  time  over  the  centralized  computation  time.  As  a  side  result, 
we  show  that  the  ACC  algorithm  may  distributively  compute  the  solution  of  any  convex  program,  and  that  it  is 
particularly  effective  when  the  dimension  of  decision  variable  is  small  compared  with  the  number  of  constraints. 
Second,  for  the  special  case  when  the  constraints  of  the  RCP  are  convex  in  the  uncertain  parameters,  we  develop 
the  vertex  constraints  consensus  (VCC)  algorithm.  We  prove  that  the  VCC  algorithm  converges  to  the  global 
solution  in  a  number  of  iterations  equal  to  the  diameter  of  the  communication  graph.  Moreover,  we  devise  a 
quantized  vertex  constraints  consensus  (qVCC)  algorithm  in  which  each  node  has  a  bounded  communication 
bandwidth.  We  provide  theoretical  bounds  on  the  number  of  the  iterations  required  for  qVCC  algorithm  to 
converge.  Last,  we  use  the  proposed  algorithms  for  distributed  classification,  distributed  estimation,  and  parallel 
model  predictive  control. 
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